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Abstract 

We present a study of S" — >■ D*^7r~and — > D*~i("^decays based on a sample of 657 x 10^ BB events collected with the 
Belle detector at the KEKB asymmetric-energy e'^e~ collider. We measure the branching fractions to be B{B^ -> D*"^7r~) = 
(1.75±0.34 (stat)±0.17 (syst)±0.11 (S)) x 10"^ and B[B° D*' K+) = (2.02±0.33 (stat)±0.18 (syst)±0.13 (B)) x 10"^ 
with significances of 6.1 and 8.0 standard deviations, respectively. The first uncertainty is statistical, the second is due to the 
experimental systematics, and the third is from uncertainties in the Dt decay branching fractions. From our measurements, 
we obtain the most precise determination so far of Rb'-k, where Rd'-r is the ratio between amplitudes of the doubly Cabibbo- 
suppressed decay B" — ^ D*^tv~ and the Cabibbo favored i?" — > D*~-k'^ decay. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 



In the standard model (SM), CP violation arises nat- 
urally when the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix is introduced into the weak-interaction 
Lagrangian A precise measurement of the CKM pa- 
rameters is crucial for understanding CP violation in the 
SM. In particular, the time-dependent CP analysis of 
the B^{B'^) D*^TT^ system provides a theoretically 
clean measurement of the product Rd-tt sin(2(/)i +(t)z) 13 j 
where (jii ^'Hd (j)^ are interior angles of the unitarity 
triangle and Rd'tt is the ratio of the magnitudes of 
the doubly Cabibbo-suppressed decay (DCSD) ampli- 
tude (Fig. ^h)) to the Cabibbo-favored decay (CFD) 
amplitude (Fig. [TJa)). Measuring the DCSD amplitude 
is not possible with the current data, due to the over- 
whelming background from D*^tt^ and hence, 
it is not possible to extract the angle 03 from this study 
alone unless an independent measurement of Rd'tt is pro- 
vided externally. The mode — >■ D*+7r'' may be used 
to estimate the size of DCSD, since B'^ — >■ D*+7r~ and 
5+ — Z?*+7r° are related by isospin symmetry How- 
ever, the 5+ — > L'*+7r° branching fraction is small and 
so far only an upper limit has been obtained Q. Unlike 
the i?° — > _D*T7r^ system, i?" _D*+7r^, which is pre- 
dominantly a spectator process with a b ^ u transition 
(Fig. [ijc)), does not have contributions from B^ decays 
to the same final state and can provide clean experimen- 
tal access to Rd*tt- Assuming SU(3) flavor symmetry 



between D* and D*, Rwtt is given by 



Rd* 



tan 0, 




B{B" D*-TT+y 



(1) 



where 9c is the Cabibbo angle, fu* and fw are the me- 
son form factors, and the i3's stand for the corresponding 
branching fractions. The D*'^it~ process, in ad- 

dition, does not have a penguin loop contribution and 
hence can in principle be used to determine \Vub\ [l]- 

In contrast to the B^ D*^tt^ decays shown in 
Figs. (Ud) and (e), the B^ £)*+7r~ decay does not 
have a contribution from the M^-exchange amplitude, as 
the quark-antiquark pair with the same flavor, required 
for such a diagram, is absent from the flnal state. We 
assume the VF-exchange contributions in B^ D*^tt^ 
to be negligible, in making the correspondence between 
L'*+7r^ and D*^it^ in the Rd*tt calculation. The size 
of the W^-exchange diagram can be estimated from the 
— )• Dl~K'^ decay, which proceeds only via W- 
exchange (Fig. [Iff)). The B" D^^ K+ branching frac- 
tion was expected to be enhanced due to rescattering 
effects m. However, a recent theoretical study based on 
measurements of related processes indicates the absence 
of such an enhancement [6| . 

While B° D+n- and B° DjK+ decays have 
been observed previously by Belle Q and BaBar , the 



observations of the modes B' 







D 



and B' 
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(c) (f) 

FIG. 1: Feynman diagrams for (a) Cabibbo-favored decay 
B° — >■ D*~7r+, (b) doubly Cabibbo-suppressed decay B° — >■ 
D'"^TT~ , (c) SU(3) flavor symmetric Ds'^-k~; color sup- 

pressed W-exchange contributions (d) to B" — >■ _D*~7r+, (e) 
to B° D'+Tv" and (f) to the decay B° D*- K+ . 

Dl~K^ have been reported by BaBar which 
measured B{B" D*+tt-) = (2.6tg| ± 0.3) x IQ-^ 
and B{B" ^ Dl~ K+) = (2.4 ± 0.4 ± 0.2) x 10"^ In 
this paper, we report an improved measurement of the 
branching fractions for the decays — >■ D*'^'k~ and 
B^ D*~K^ with a data sample consisting of 
657 X 10^ BB pairs, collected with the Belle detector 
at the KEKB asymmetric-energy e+e" collider [T2| . 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector, a 50- 
layer central drift chamber (CDC), an array of aerogel 
threshold Chercnkov counters (ACC), a barrel- like ar- 
rangement of time-of-flight scintillation counters (TOF) , 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux- return 
yoke located outside the solenoid is instrumented to de- 
tect Kl mesons and to identify muons. The detector is 
described in detail elsewhere [l3|. Two different inner 
detector configurations were used. For the first sample 
of 152 X 10^ BB pairs, a 2.0 cm radius beam-pipe and 
a 3-layer silicon vertex detector were used; for the latter 
505 X 10^ BB pairs, a 1.5 cm radius beam-pipe with a 
4-layer silicon vertex detector and a small-cell inner drift 
chamber were used [l^l . 

The signal is reconstructed in three modes: (pn^ 
with (j) K+K~, K*{SmfK+ with 7?* (892)° ^ 
K~t:+, and K^K^ with tt+tt". Charged tracks 

arc selected with requirements based on the impact pa- 
rameter relative to the interaction point (IP). The devia- 
tions from the IP are required to be within ±4 cm along 



the 2;-axis (the direction opposite to the positron beam) 
and within 0.2 cm in the x — y plane. We also require the 
transverse momentum of the tracks to be greater than 0.1 
GeV/c in order to reduce low momentum combinatorial 
background. 

For charged particle identification (PID), we combine 
the information from the specific ionization (dE/dx) in 
the CDC with measurements from the TOF and ACC. 
At large momenta {p > 2.5 GcV/c) only the ACC mea- 
surement and dE/dx arc used. Wc assign likelihood 
values Ck (^tt) for the kaon (pion) hypothesis to each 
charged track. Tracks are identified based on the ratio 
T^k/tt = C,kI{^k + -C^), which peaks at one for real 
kaons and at zero for real pions. For the prompt kaon 
(pion) track, we require 'R-Kj-w > 0.6(< 0.6), for which the 
identification efficiency is 85% (92%) with a pion (kaon) 
fake-rate of 8% (15%). Due to the low background level 
for modes with a <p meson, less restrictive PID cuts of 
T^k/tt > 0.2 are applied to the (j> daughter tracks. 

The 4> (A'*(892)°) mesons are required to have an 
invariant mass within ±14 MeV/c^ (±75 MeV/c^) of 
the nominal (j) {K*{892)^) mass [31 . We reconstruct 
Kg candidates from tt+tt^ pairs, requiring the invariant 
mass to be within ±10 MeV/c^ (~ ±3cr) of the nomi- 
nal Kg mass. The Kg candidate is further required to 
pass a momentum-dependent selection criteria based on 
its vertex topology, the flight length in the r — plane, 
and the daughter tt* momentum distribution |16| . The 
Df candidate mass window for (^7r+, X*(892)"A'+, and 
K^K+ modes is ±13 MeV/c^ ±15 MeV/c^, and ±17 
MeV/c^, respectively; these ranges correspond to approx- 
imately 3ct in resolution around the mass. To reduce 
combinatorial background, we use a more stringent PID 
requirement for the kaon accompanying the A"* (892)°, 
TZk/tt > 0.8. The candidate is constrained kinemat- 
ically to have a mass equal to the nominal value (Tsj . 

The D*^ mesons arc reconstructed by combining the 
Df candidates with a photon. The photons are recon- 
structed from energy depositions in the ECL and are re- 
quired to have energies greater than 60 MeV (100 MeV) 
in the barrel (endcap) region covering the polar angle 
32° < 6* < 128° (17° < 6* < 32° (forward endcap) and 
128° < e < 150° (backward endcap)). The can- 
didate is required to have AM = Mj^+^ — Mjj+ be- 
tween 128 MeV/c^ and 162 MeV/c^, where M^+^ and 
Mjj+ are the invariant masses of the Dfj system and 
the Df candidate, respectively. To reduce the combina- 
torial background due to low energy photons, we require 
that cos 61^.+ > -0.6 (-0.7) for B° D*+tt- {B° 
D*^K^), where 9jy+ is defined as the angle between the 
flight direction of the photon and the direction opposite 
to the B° flight in the D*~^ rest frame. We then per- 
form a mass-constrained flt to the candidate. This 
improves the momentum resolution by 25%. 

The candidates, reconstructed by combining a D*^ 
candidate with an oppositely charged pion/kaon track, 
arc identified by the energy difference, AE = J^i ~ 
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Ehcam and the beam-energy constrained mass, Mbc = 
V^boam - (ZiPi)^; where Sboam IS the beam energy in 
the T(45') center-of-mass (CM) frame and pi and Ei are 
the momentum and energy of the ith daughter of the 
in the CM frame. We retain 5° candidates with AE 
within ±0.2 GcV and Afbc between 5.2 GeV/c^ and 5.3 
GeV/c^ for further analysis. 

The dominant background comes from the e'^e" — >■ qq 
(q = u,d, s and c quarks) continuum process. To sup- 
press this background, we use the event topology in 
the CM frame to distinguish more spherical BB events 
from the jet-like continuum events. A likelihood func- 
tion TZ = £sig/('Csig + -Cbkg) is prepared by combining 
a Fisher discriminant, based on a set of modified Fox- 
Wolfram moments [U with cos 9b, where 9b is the 
polar angle of the B^ meson flight direction in the CM 
frame. The angle 9b follows a sin^^s distribution for 
BB events, while the continuum distribution is flat. The 
selection criteria for TZ arc determined by maximizing a 
figure-of- merit, S/y/S + B, where S and B arc the num- 
ber of signal and background events determined from 
large Monte Carlo (MC) samples [l^, with statistics cor- 
responding to about 100 (5) times data for signal (back- 
ground) MC. The signal yield, S is obtained assuming 
the latest branching fraction measurements [l^. In the 
case of B° D*+it- {B° D*-K+), we require U to 
be greater than 0.45 (0.45) for the (jiir mode, 0.50 (0.60) 
for the K*{S%2fK mode, and 0.40 (0.40) for the K^K 
mode. For the 07r mode, this requirement suppresses 80% 
of the continuum background, while retaining 85% of the 
signal. 

About 15% of events have more than one B^ signal can- 
didate. For these events we choose the candidate with the 
Mbc value closest to the nominal B'^ mass. This proce- 
dure selects the correct B^ candidate in about 92% of the 
cases. Only events with Mbc between 5.27 GeV/c^ and 
5.29 GeV/c^ are considered for further analysis, while 
the signal is extracted by performing a fit to the /S.E 
distribution. We define the fit region to be \A.E\ < 0.2 
GeV. 

A MC sample of BB events is used to determine pos- 
sible backgrounds that can enter the AE fit region. In 
both signal modes (i?° — > D*h, where h is a charged K 
or tt), about 45% of the background comes from decays 
involving a K^tt^tt^ or a Z)+ ^s""^ 

decay, which form a fake Df when the 7r+ from the Z?^ 
is misidcntificd as a K~^. However, these events do not 
peak and are distributed over the entire fit region due to 
the addition of a random photon. 

On the other hand, some rare B decays to final states 
that contain a correctly reconstructed Ds produce 
non-negligible peaking structures in the AE fit region: 
B° D+TT- {B° events populate the re- 

gion around 150 MeV due to the addition of an extra 
photon, B° D*+p- {B+ D*-K+TT+) events pop- 
ulate the region around —150 MeV, since a 7r° (tt"*") is 
not reconstructed, while the events from i?° Dfp^ 
(J5+ DjK+'K+) are distributed around -50 MeV, 



with the extra photon compensating the lost pion. These 
backgrounds are represented by PDFs with fixed yields. 
MC samples are used to determine the PDF parameters 
as well as efficiencies. The yields are calculated assuming 
the most recent known values for their branching frac- 
tions [l3, [l^. Apart from the backgrounds discussed 
above, the two signal modes cross-feed each other. 
We use signal MCs to determine the PDFs for the cross- 
feeds. AM sidebands in the data are used to verify the 
consistency of the MC background predictions with the 
data. 

The signal PDF is the sum of a Crystal Ball line- 
shape [20| and a broad Gaussian and is parametrized 
using signal MC samples generated in each Df mode. 
The signal as well as the peaking background PDFs 
are subsequently corrected for possible differences in the 
parameter vaues between MC and real data, using a 
B^ — )• Dl~^D~ data control sample. The combinatorial 
backgrounds in each mode are accounted for by adding 
linear functions. 

We determine the branching fractions from a simulta- 
neous unbinned extended maximum likelihood fit to the 
AE distributions of the three Df decay modes for each 
signal mode. To account for the cross-feeds between the 
signal modes due to the misidentification of the prompt 
track, the two B signal modes are fitted simultaneously, 
with the — )- D*~K~^ signal yield in the correctly re- 
constructed sample determining the normalization of the 
cross-feed in the B'^ — >■ 7r~flt region, and vice versa. 
The fit has 14 free parameters: the branching fractions 
of the signal modes (2), and the yields and slopes of 
the first-order polynomials representing the combinato- 
rial background in each of the three modes (12). Fig- 
ure [H shows results of the simultaneous fit for the three 
Df modes in both signal modes. We summarize the re- 
su lts of the fits in Table [J The significance is defined as 
■\/— 2 ln(£o/'Cmax), where £max (Co) are the likelihoods 
for the best fit and with the signal branching fraction 
fixed to zero. 

The major source of systematic uncertainty in the 
branching fraction measurement of B'^ — )- D*~^7t~ {B'^ — >■ 
Dl~K'^) is the uncertainty in the branching fractions 
of the D+ decays, which amount to 5.9% (6.2%). The 
uncertainties in the branching fractions of the peaking 
background modes contribute an additional error of 1.5% 
(1.9%). The systematic uncertainty in the tracking effi- 
ciency is estimated to be about 1.0% per track. The 
uncertainty in the PID efficiency is about 2.4% (2.1%). 
Photon detection efficiency has an uncertainty of 7.0%, 
while Kg detection efficiency adds 1.1% uncertainty in 
the result. The efficiency of the TZ requirement, used to 
suppress the continuum background, introduces an un- 
certainty of 0.6% (0.5%) in the branching fraction. The 
limited size of the MC samples used to determine effi- 
ciencies and cross-feed fractions introduces an error of 
1.4% (1.6%). The uncertainty in the determination of 
the signal PDF shape is about 3.4% (1.5%). Estimation 
of possible bias in the fit results in another 0.9% (0.3%) 
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TABLE I: Efficiency (e), yield (A'sig), branching fraction (B) and statistical significance not including systematic errors (5) 
from the fits to the data obtained individually in the three Df modes as well as from the simultaneous fit. The second error 
on the B's is due to the uncertainties in Df decay branching fractions. The individual fit results are consistent with each other 
and also with the simultaneous fit. Note that the efficiencies in the K%K mode include the K% — > tt^tt" decay branching 
fraction. 
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3.2 


5° 








7.9 
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1.78 ±0.92 ±0.11 
1.75 ±0.34 ±0.11 
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23 ±7 
14 ±5 
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2.02 ±0.33 ±0.13 
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TABLE IL Contributions to the systematic uncertainty 
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FIG. 2: The simultaneous fit in the B° -i> D*+7r" ((a) cf>Tv, (b) 
K*°K and (c) K'^K mode) and the B° D*' K+ ((d)-(f)) 
signal modes. Signal peaks are shown by the solid curves while 
the solid-filled curves represent the cross-feed contributions 
from the other signal modes. The long-dashed curves cor- 
respond to the contribution from the B — i- Dg-K {B — ^ DgK) 
and the dot-dashed curves to that from 



D. 



(*)+, 



{B+ Dr>'K+n+). The dotted curves correspond to the 
combinatorial background. 



Source 



Contribution(%) 
D*+Tr- D*+K- 



Ds branching fraction uncertainties 



signal 


5.9 


6.2 


peaking background 


1.5 


1.9 


Total (S) 


6.1 


6.5 


Tracking efficiency 


4.0 


4.0 


Photon detection efficiency 


7.0 


7.0 


Particle identification efficiency 


2.4 


2.1 


K% efficiency 


1.1 


1.1 


cn 


0.6 


0.5 


Nbb 


1.4 


1.4 


MC statistics 


1.4 


1.6 


PDF shape 


3.4 


1.5 


Fit bias 


0.9 


0.3 


Total (other) 


9.4 


8.8 



0.17 (syst) ± 0.11 (B)) X 10"^ and B{B^* D*- K+) = 
(2.02 ± 0.33 (stat) ± 0.18 (syst) ± 0.13 (6)) x 10"^ with 
significances of G.lcr and 8. Oct, respectively, where the 
systematic uncertainties on the signal yield as well as 
the statistical uncertainties are included in the signif- 
icance evaluation. Though consistent with the previ- 
ous measurements [To| , wc observe slightly lower branch- 
ing fractions. Using the observed value for the i?" ^■ 
£'*+7r~ branching fraction, the latest values for B{B^ — > 
D*-t:+) ^ (2.76 ± 0.13) x 10"^^ tan6ic = 0.2314 ± 
0.0021 [l^, and the theoretical estimate of the ratio 
fjj+/fD+ = (1.164 ± 0.006 (stat) ± 0.020 (syst)) [H, 
we obtain, 



uncertainty. 

Table |TT] summarizes the systematic uncertainties in- 
volved. The overall systematic error is obtained by 
adding the above contributions in quadrature. 

Wc obtain B{B° ^ D*+it^) = (1.75 ± 0.34 (stat) ± 



Rij,^ (1.58 ± 0.15(stat) ± O.lO(syst) ± 0.03(th))%, 

where the first error is statistical, the second corresponds 
to the experimental systematic uncertainty and the third 
accounts for the theoretical uncertainty in the /^+ //£>+ 
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estimation. We have assumed that the ratio /dJ/d 
is equal to the ratio of vector meson decay constants, 
fo' / fo' ■ The quenched QCD approximation 12^ as well 
as the heavy quark effective theory predictions j23| point 
toward an uncertainty of about 1% due to this assump- 
tion, which is included in our estimation of Rd-tt- The 
value we obtain for Rd-'tt, though consistent with the 
theoretical expectation of 2%, is slightly smaller than the 
previous estimate [lOj . 

The observed value for the i?° — >■ D*~K^ branch- 
ing fraction is two orders of magnitude lower than that 
for the Cabibbo-favored decay — >■ D*~tt~^. This can 
be understood purely in terms of the exchange ampli- 
tude and there is no evidence for enhancement due to 
rescattering effects, which would lead to comparable am- 
plitudes for the two processes Q. From this same com- 
parison, we find no evidence for large M^-exchange contri- 
butions to -8° — >■ D*^TT^; such contributions arc assumed 
to be small, in the determination of Rd'tv (Eq. [T]). 

In conclusion, we report the most precise measurement 
of the B° D*+-K- and B° -J> D*-K+decscy branching 
fractions. This improves the precision with which the 
parameter Ro'-n can be estimated, and thus the prospect 
of determining <f)^ from measurements of CP violating 
effects in the D*^tt^ system. 
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